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ABSTRACT 


This  paper  presents  the  results  of  an  investigation 
of  the  flutter  characteristics  of  a  fully  submerged 
subcavitating  hydrofoil.  The  foil  suffered  a  cata¬ 
strophic  flutter  failure  at  a  speed  of  approximately 
48  knots.  Comparisons  are  made  with  various  pre¬ 
diction  methods,  neither  the  two-dimensional  nor 
the  three-dimensional  theoretical  calculations 
yielding  results  of  the  correct  order  of  magnitude. 
Calculations  employing  measured  three-dimensional 
oscillatory  coefficients  give  a  highly  conservative 
prediction  for  flutter  speed. 


INTRODUCTION 


The  possible  occurrence  of  classical  bending-torsion 

flutter  of  a  fully  submerged  subcavitating  hydrofoil  has  been 

questioned  repeatedly  over  the  past  several  years.  In  the 

absence  of  definitive  and  authoritative  evidence  of  such  a 

flutter  occurrence,  or  related  experimental  data,  opinions 

have  ranged  from  the  impossibility  of  flutter  to  speculation 

as  to  the  inapplicability  of  conventional  aerodynamic  flutter 

theory.'*'  As  part  of  a  fundamental  experimental  investigation 

into  these  questions,  the  authors  have  conducted  a  series  of 

tests  with  special  hydrofoil  models  capable  of  providing 

directly  the  three-dimensional  coefficients  of  oscillatory  lift 

and  moment  due  to  cantilever  wing  bending  and  torsion,  in  water. 

This  data  covers  a  range  of  low  values  of  reduced  velocity  and 

2,3 

has  been  presented  in  detail  elsewhere.  A  bending-torsion 
flutter  model  was  also  constructed  and  tested  as  part  of  this 
research  program,  and  it  is  the  presentation  of  data  regarding 
this  flutter  model  to  which  the  present  paper  is  directed. 


EXPERIMENTAL  PROGRAM 


Description  of  Flutter  Model 


The  hydrofoil  flutter  model  (Fig.  1)  was  deliberately 


designed  to  be  flutter  prone,  in  order  to  provide  a  definite  experi¬ 
mental  point  for  correlation  with  various  analytical  predictions. 

Vr>  . 

A  mass  ratio  (  —  —  ■g 4  -  )  of  unity  was  arbitrarily  selected  for 

the  model  as  being  small  enough  to  demonstrate  behavior  in  the 
range  of  interest  to  hydrofoil  applications,  but  large  enough  to 
insure  that  flutter  would  occur.  This  mass  ratio  is  also  about 
as  large  a  value  as  can  be  practically  achieved  with  a  12%  thick 
foil  section  of  solid  lead.  An  NACA  16-012  airfoil  section  was 
chosen  in  order  to  have  the  center  of  gravity  near  midchord. 

The  elastic  axis  was  located  well  ahead  of  the  c.  g.  ,  and  the 
bending  and  torsional  stiffnesses  were  held  low  in  order  to  ob¬ 
tain  a  low  flutter  speed.  Further,  in  order  to  simplify  and 
enhance- the  accuracy  of  the  analytical  procedures,  the  model 
parameters  were  made  uniform  spanwise  and  the  spar,  being 
the  only  continuous  spanwise  member,  provided  a  well  defined 
elastic  axis. 

The  remaining  model  design  parameters  were  also  chosen 
on  the  basis  of  values  that  would  result  in  the  lowest  practicable 
flutter  speed,  as  given  in  Table  I,  The  actual  values  of  these 
various  parameters,  as  determined  by  measurements  from  the 
model,  are  also  shown.  While  the  uncoupled  bending  frequency 
was  found  to  be  close  to  the  design  value,  the  torsional  frequency 
was  approximately  twice  as  large  as  the  design  value,  primarily 


as  the  result  of  failing  to  account  properly  for  the  stiffening  effects 
of  the  spar  attachments. 

The  model  was  composed  of  isolated  spanwise  strips  con¬ 
nected  by  a  single  H-beam  spar.  The  spar  was  machined  from 
4340  steel  and  was  designed  to  have  a  high  chordwise  bending 
stiffness  as  compared  with  the  lateral  bending  and  torsional  stiff¬ 
ness.  The  foil  contour  was  formed  by  cast  lead  sections  with 
embedded  steel  webs  to  provide  additional  rigidity.  The  webs 
were  first  welded  to  steel  blocks,  which  then  were  allowed  to 
project  from  the  lead  castings  to  provide  for  mounting  to  the 
spar  caps.  A  typical  section  is  hown  in  Figure  2.  Each  lead 
section  was  approximately  1.  90  inches  wide,  separated  from  the 
adjacent  section  by  a  gap  of  approximately  0.10  inch  and  sealed 
with  Silastic  RTV  731  approximately  0.10  inch  thick.  The  foil 
contour  was  filled  out  over  the  spar  caps  with  a  lead-filled  plastic 
poured  in  place  and  trimmed  to  contour.  This  material  was  of 
such  composition  that  it  would  hold  its  shape,  but  without  contri¬ 
buting  significantly  to  the  spar  stiffness.  A  stainless  steel  plate 
was  employed  to  close  the  tip  of  the  model. 

Model  Suspension  and  Excitation  System 

The  spar  entended  past  the  model  root  and  fittings  were 

attached  for  mating  with  the  same  support  fairing  used  for  the  load- 
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measuring  models  employed  earlier.  The  root  end  was  therefore 


held  rigidly  so  that  the  model  was  cantilevered  vertically  from  the 
towing  carriage  (Fig.  1). 

An  internal  system  for  dynamic  excitation  of  the  model  was 
desired  in  order  that  no  disturbance  to  the  fluid  flow  would  result. 
This  system  consisted  simply  of  two  7 / 16  in.  diameter  rods  passing 
through  the  spar  on  each  side  of  the  web  (Fig.  2)  and  anchored  to 
the  model  tip,  but  free  floating  along  the  rest  of  the  span  and  ex¬ 
tending  through  the  model  root  and  up  to  the  upper  side  of  the 
carriage  frame.  Crank  arms  connected  with  links  were  attached 
to  the  upper  ends  of  the  rods  and  a  cable  attached  to  the  end  of 
these  cranks  extended  to  the  sides  of  the  towing  carriage.  When 
the  cable  was  suddenly  pulled,  the  torque  applied  to  the  rods 
twisted  the  model  tip  slightly.  The  rods  were  also  designed  to 
serve  as  components  of  a  brake  system  to  prevent  structural 
damage  to  the  model  should  a  violent  mode  of  flutter  be  encountered. 
For  this  purpose,  the  cranks  at  the  top  ends  of  the  rods  could  be 
clamped  between  two  manually  operated  levers  to  increase  the 
overall  torsional  stiffness  of  the  model.* 


*  Because  the  torsional  stiffness  of  the  model  was  actually  much 
higher  than  originally  anticipated,  the  added  stiffness  provided 
by  the  clamped  rods  was  not  significant  and  thus  the  braking 
system  was  rendered  ineffective. 


Instrumentation  and  Test  Procedures 


Instrumentation  consisted  of  two  strain  gage  bridges  on 
the  spar  root,  one  being  sensitive  to  spar  bending  and  one  to 
spar  torsion.  Signals  from  these  bridges  were  recorded  on  an 
oscillograph  and,  additionally,  the  torsion  signal  was  also  dis¬ 
played  on  an  oscilloscope.  A  high-speed  motion  picture  camera 
was  stationed  at  a  single  point  along  the  length  of  the  towing  basin. 
All  tests  were  conducted  on  the  high  speed  carriage  at  the  David 
Taylor  Model  Basin. 

The  towing  carriage  speed  was  varied  in  increments  of 
approximately  5  knots,  from  an  initial  speed  of  10  knots.  As  soon 
as  a  given  speed  was  established,  the  recording  oscillograph  was 
started  and  the  exciting  cable  pulled,  so  that  logarithmic  decrements 
of  the  signal  decay  from  the  two  strain  gage  bridges  were  thereby 
recorded.  The  decay  of  the  torsion  signal  was  simultaneously 
observed  on  the  oscilloscope  (Polaroid  photo)  and  the  logarithmic 
decrement  obtained  immediately.  Thus,  it  was  possible  to  maintain 
a  continuing  surveillance  of  the  system  damping  as  carriage  speed 
increased.  Some  of  the  records  taken  at  the  higher  carriage  speeds 
were  repeated. 

Catastrophic  Failure 


At  a  carriage  speed  of  approximately  45  knots,  the  measured 


logarithmic  decrement  showed  a  significant  decrease  in  damping 
and  therefore  the  next  speed  increment  was  selected  to  be  only 
2-1/2  knots.  However,  the  actual  carriage  speed  was  48.1  knots, 
resulting  in  a  strong  flutter  condition  growing  from  the  initial 
disturbance  applied  to  the  model  and  leading  to  catastrophic 
failure  of  the  model  within  a  few  seconds.  Structural  failure 
occurred  at  the  spar  root,  as  shown  in  Figure  3.  The  model 
then  slid  off  of  the  exciter  rods  (which  were  bent  to  an  angle  of 
about  90°),  rose  completely  out  of  the  water,  and  impaled  itself 
in  the  rear  wall  (1/8  in.  aluminum  sheet)  of  the  test,  bay  of  the 
carriage  (Fig.  4). 

The  flutter  frequency  was  approximately  17.  5  cps,  resulting 
in  a  value  of  reduced  velocity  at  flutter  of  (V/bw)j,  =  1.  48.  By  a 
rather  fortuitous  circumstance,  the  oscillatory  motion  of  the  foil 
had  just  attained  rather  large  amplitudes  as  it  passed  the  point  at 
which  the  motion  picture  camera  was  stationed  so  that  a  short 
film  record  was  obtained.  The  film  shows  some  cavitation  arising 
from  these  large  amplitude  motions,  just  prior  to  the  structural 
failure. 


ANALYSIS 

Three  flutter  analyses  have  been  performed  in  order  to 


compare  various  prediction  methods  with  the  measured  flutter 
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speed  of  this  model.  The  various  analyses  are  all  based  on  a 

ten-station  strip  theory  Rayleigh-type  calculation^,  and  differ 

from  each  other  only  in  the  values  employed  for  the  oscillatory 

lift  and  moment  coefficients. 

One  analysis  was  made  employing  conventional  two- 

4 

dimensional  theoretical  coefficients  and  one  was  made  employ- 

5 

ing  modified  three-dimensional  theoretical  coefficients;  the 
third  analysis  employed  the  measured  three-dimensional 
coefficients  .  The  measured  three-dimensional  coefficients 
were  originally  obtained  in  terms  of  moment  about  the  31%  chord 
position.  The  more  conventional  quarter -chord  coefficients  are 
then  obtained  from  the  relations 


Lh  = 

L  hh 

- 

L  hit 

Mh  = 

M  oth  +  mb/b  Lhh 

M*  - 

+  mt/b  Lhh 

where  m^/b  =  0. 1131  for  the  bending  model  and  mj/b  =  0.  1108  for 
the  torsion  model.  The  elastic  axis  was  located  at  the  quarter- 
chord  position  in  both  models.  Values  of  the  quarter-chord 
coefficients  for  the  ten  spanwise  stations  and  nine  nominal  values 
of  reduced  velocity  are  given  in  Tables  II-X,  as  interpolated 
from  the  actual  experimental  values. 


The  results  of  the  various  flutter  analyses,  together  with 
the  relevant  experimental  data*  discussed  earlier  in  this  paper, 
are  presented  in  Figure  5  in  terms  of  the  exponential  decay 
factor  (&)  versus  speed  (V).  Points  at  which  these  curves 
intercept  the  speed  axis  are  defined  as  critical  flutter  speeds. 

Another  comparison  between  the  various  analyses  and  the 
experimental  data  is  shown  in  Figure  6  in  terms  of  the  coupled 
torsional  frequency  versus  speed.  The  calculated  values  employ¬ 
ing  theoretical  three-dimensional  coefficients  appear  to  give 
a  trend  most  closely  approximating  that  of  the  experimental  data. 

DISCUSSION 

The  results  given  in  Figure  5  are  striking  in  their  dif¬ 
ferences.  The  two-dimensional  theory  predicts  no  critical  speed 
at  all  while  the  three-dimensional  theory  leads  to  a  critical  speed 
that  is  unconservative  by  more  than  100%,  The  calculation  em¬ 
ploying  measured  three-dimensional  coefficients,  on  the  other 
hand,  predicts  a  flutter  speed  that  is  conservative  by  a  substantial 
amount.  Obviously,  none  of  these  can  be  considered  a  satisfactory 

*  The  multiple  experimental  points  at  certain  values  of  speed  show 
the  scatter  of  values  of  logarithmic  decrement  obtained  for  dif¬ 
ferent  test  runs  at  that  speed. 


prediction  of  critical  flutter  speed. 

Based  on  previous  aeronautical  experience  with  similar 
techniques,  ^  it  would  normally  be  expected  that  the  calculation  of 
flutter  speed  employing  measured  three-dimensional  coefficients 
would  yield  a  value  in  very  close  agreement  with  the  measured 
flutter  speed.  As  was  pointed  out  in  Reference  3,  however,  the 
measured  coefficients  employed  here  are  known  to  be  somewhat 
low  in  magnitude,  particularly  the  moments,  due  to  attenuation 
of  response  of  the  measuring  system  to  loads  aft  of  midchord, 
with  an  unknown  effect  on  phase  angles.  This,  coupled  with  some 
uncertainties  in  data  reduction  due  to  impure  modes  and  high  noise 
levels  lends  some  doubt  as  to  the  ability  to  predict  accurate  flutter 
speeds  with  these  measured  coefficients. 

Another  factor  which  may  warrant  consideration  resides 
in  the  observation  of  some  partial  cavitation  occurring  on  the  model 
during  the  large  amplitude  oscillations  just  preceding  structural 
failure.  There  is,  however,  no  evidence  that  this  cavitation 
occurred  before  large  amplitudes  were  built  up. 

Perhaps  one  more  observation  may  be  pertinent  to  this 
discussion.  During  the  course  of  these  calculations,  some  small 
changes  in  the  numerical  values  of  the  various  coefficients  were 


made  with  very  large  effects  on  the  value  of  the  critical  flutter 
speed..*  At  least,  one  should  conclude  that,  for  the  model 
parameters  involved  here,  the  critical  flutter  speed  is  extremely 
sensitive  to  the  values  of  the  oscillatory  coefficients. 

The  results  of  analytical  predictions  of  flutter  speeds  of 
fully  submerged  subcavitating  hydrofoils  should  be  viewed  with 
skepticism,  and  substantiating  evidence  from  model  tests  should 
be  obtained  whenever  possible. 


These  studies,  aimed  at  developing  a  semi-empirical  flutter 
prediction  technique  more  generally  valid  than  the  conven¬ 
tional  calculations  reported  here,  are  continuing. 
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TABLE  I.  FLUTTER  MODEL  DESIGN  PARAMETERS 


Model  Design  Parameters 


Design  Value 


Actual  Value 


Aspect  ratio, 

Semichord,  b 

Mass  ratio, yit 

Elastic  axis  location,  a 

Center  of  gravity  location,  xM 

Radius  of  gyration,  r* 

Bending  stiffness,  El 
Torsional  stiffness,  GJ 
Frequency  ratio,  ^fuJx 
Torsion  frequency,  u3« 

Total  weight 


5.  00 


0.  50  ft. 
1.  00 
-  0.  50 


0.  50 
0.  50 

3.  08  x  106 
0.  311  xlO6 
1.  00 

10. 0  cps* 


lb-in? 

lb-in 


5.00 


0.  50  ft. 
0.  99 
-  0.  50 


0.  524 
0.  512 

3.40  x  106  ib-in2 
0.973  x  106  ib-in2 
0.  490 
20. 5  cps* 

121.  2  lb. 


*  Uncoupled  natural  frequencies  in  air 


TABLE  II.  OSCILLATORY  HYDRODYNAMIC  COEFFICIENTS  FOR  l/k=  0 


Span  Mag.  Phase  Mag.  Phase  Mag.  Phase  Mag.  Phase 


Mag. 

Phase 

.  561 

0 

.  376 

0 

.  315 

0 

.  268 

0 

.  232 

0 

.  198 

0 

.  173 

0 

.  150 

0 

.  128 

0 

0 

0 

TABLE  III.  OSCILLATORY  HYDRODYNAMIC  COEFFICIENTS  FOR  l/k=.  500 


Phase 


Phase  Mag.  Phase  Mag.  Phase 


8.  5° 
13.  3° 
17.  3° 
20.  9° 
23.  5° 
25.8° 
28.  1 0 
30.  2“ 
31.3° 
32.  0° 


TABLE  IV.  OSCILLATORY  HYDRODYNAMIC  COEFFICIENTS  FOR  l/k=  .700 


% 

Lh 

M  h 

M* 

Span 

Mag. 

Phase 

Mag. 

Phase 

Mag. 

Phase 

Mag. 

Phase 

10 

1.  84 

1. 467 

251° 

.  567 

10.  1° 

.  715 

20 

1. 28 

311.5° 

1.016 

256° 

.  462 

17.  8° 

.  516 

30 

.  905 

316° 

.  995 

260° 

.  393 

24.  3° 

.  468 

311.5° 

40 

.  729 

.949 

263.  5° 

.  341 

IV 

OO 

o 

© 

.  426 

50 

.  668 

325  ° 

.  886 

267° 

.  306 

31.1° 

.  393 

315° 

60 

.  667 

329° 

.  824 

269° 

.  283 

32.  1“ 

316° 

70 

.  658 

331° 

271° 

.  266 

33.  7° 

.  334 

316° 

80 

.  628 

334° 

.  704 

272. 5° 

.  247 

35.  5° 

.  314 

317° 

90 

.  531 

335° 

.  626 

.  199 

36.4° 

.  295 

317° 

100 

0 

336° 

273° 

0 

37.0° 

0 

317° 

TABLE  V.  OSCILLATORY  HYDRODYNAMIC  COEFFICIENTS  FOR  l/k=  .833 


% 

h 

L* 

Mh 

M* 

A 

Span 

Mag. 

Phase 

■ 

Phase 

Mag. 

Phase 

10 

2.210 

300° 

2.  00 

240° 

.  589 

11.0° 

.  797 

299“ 

mm 

1.515 

304° 

1. 49 

245.  5° 

.  467 

18.  6° 

.  590 

303° 

■9 

1. 065 

308. 5° 

1.  38 

250° 

.  387 

25.  4° 

.  537 

306“ 

40 

.  850 

313° 

1.  27 

254.  5° 

.  333 

29.2“ 

.  496 

309“ 

50 

.  765 

318° 

1.  16 

258. 5° 

.  297 

32.  5° 

.  460 

312“ 

60 

321. 5° 

i.  06 

261  ° 

.  274 

33.  7° 

.  427 

313° 

70 

.  740 

324° 

.  975 

263. 5° 

.  258 

35.  6° 

.  397 

315° 

80 

.  695 

327° 

.  890 

265° 

.  243 

38.  1  ° 

.  373 

316“ 

90 

.  570 

329° 

265° 

.  197 

39.4° 

.  347 

317° 

100 

0 

265° 

40.  0“ 

0 

317° 
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FIGURE  2  -  DETAILS  OF  FLUTTER  MODEL  CONSTRUCTION 
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